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A search for 21 cm H i absorption in AT20G compact radio 
galaxies 



(N 

o 

U 

Or 
< 

(N 



o 

u 

Ok 
6 

OS 

m 



J. R. Allison 1 *, S. J. Curran 1 ' 2 , B. H. C. Emonts 3 , K. Gereb 4 ' 5 , E. K. Mahony 1,2 t, 
S. Reeves 1 ' 2 , E. M. Sadler 1 ' 2 , A. Tanna 6 , M. T. Whiting 3 , and M. A. Zwaan 7 

^Sydney Institute for Astronomy, School of Physics A28, University of Sydney, NSW 2006, Australia 

2 ARC Centre of Excellence for All-sky Astrophysics ( CAASTRO) 

3 CSIRO Astronomy & Space Science, P. O. Box 76, Epping NSW I 710, Australia 

4 ASTRON, The Netherlands Institute for Radio Astronomy, Posthus 2, 7990 AA, Dwingeloo, The Netherlands. 

5 Kapteyn Astronomical Institute, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands 

6 School of Physics, University of New South Wales, Sydney NSW 2052, Australia 

7 European Southern Observatory, Karl-Schwarzschild-Str. 2, Garching 85748, Germany 



ABSTRACT 

We present results from a search for 21 cm associated H I absorption in a sample of 
29 radio sources selected from the Australia Telescope 20 GHz survey. Observations 
were conducted using the Australia Telescope Compact Array Broadband Backend, 
with which we can simultaneously look for 21cm absorption in a redshift range of 
0.04 < z < 0.08, with a velocity resolution of 7 km s _1 . In preparation for future large- 
scale H I absorption surveys we test a spectral-line finding method based on Bayesian 
inference. We use this to assign significance to our detections and to determine the 
best-fitting number of spectral-line components. We find that the automated spectral- 
line search is limited by residuals in the continuum, both from the band-pass calibration 
and spectral-ripple subtraction, at spectral-line widths of Aiifwhm ^ 10 3 kms _1 . Us- 
ing this technique we detect two new absorbers and a third, previously known, yielding a 
10 per cent detection rate. Of the detections, the spectral-line profiles are consistent with 
the theory that we are seeing different orientations of the absorbing gas, in both the host 
galaxy and circumnuclear disc, with respect to our line-of-sight to the source. In order 
to spatially resolve the spectral-line components in the two new detections, and so verify 
this conclusion, we require further high-resolution 21 cm observations (~ 0.01 arcsec) 
using very long baseline interferometry. 
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1 INTRODUCTION 

Observations of the 21 cm Hi line of neutral atomic hydrogen 
in absorption against bright radio continuum sources can pro- 
vide a unique probe of the distribution and kinematics of gas 
in the innermost regions of radio galaxies. For a fixed back- 
ground continuum source brightness the detection limit for H I 
and molecular absorption lines is independent of the redshift of 
the absorber. Therefore measurements of associated H I absorp- 
tion in active galaxies are sensitive to relatively small amounts 
of neutral gas, provided the central radio continuum source is 
sufficiently strong and compact. Hi absorption-line measure- 
ments can detect both circumnuclear gas discs and large scale 
gas, including inflows and outflows related to to Active Galactic 
Nuclei (AGN) fuelling and feedback (e.g. Morganti et al. 2001; 
Holt, Tadhunter, & Morganti 2008; Morganti et al. 2011). Re- 
cent studies by Morganti and colleagues have found broad, shal- 
low absorption-line components arising from the large-scale gas 
outflows in radio galaxies (Morganti, Tadhunter, & Oosterloo 
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2005; Holt et al. 2008), which provide some of the most direct 
evidence for AGN-driven feedback in massive galaxies (Croton 
et al. 2006). 

Previous work may suggest that the detection rate of associ- 
ated H I absorption is highest in the most compact radio galaxies, 
Morganti et al. (2001) searched for H I absorption in 23 nearby 
(z < 0.22) radio galaxies from the southern 2 Jy sample (Mor- 
ganti et al. 1999), from which they obtained five detections (an 
overall detection rate of 22 per cent). They found that the Hi 
detection rate in compact radio galaxies was much higher than 
in classical FR-I and broad emission-line FR-II systems. Fur- 
thermore Vermeulen et al. (2003) detected H I absorption in 19 
(33 per cent) of a sample of 57 compact radio galaxies at z < 
0.85. Pihlstrom, Conway, & Vermeulen (2003) show that the in- 
tegrated optical depth for 21 cm Hi absorption increases with 
decreasing source size, suggesting that the probability of detec- 
tion might be highest for the more compact sources. However 
Orienti, Morganti, & Dallacasa (2006) show that this relation- 
ship might break down for the most compact (high-frequency 
peaker) sources. Recent searches for associated H I absoiption 
in nearby galaxies have also found relatively high numbers of 
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detections by targetting compact radio sources (e.g. Gupta et al. 
2006; Chandola, Sirothia, & Saikia 201 1). 

In the work presented here we use the recently-completed 
Australia Telescope 20 GHz survey (AT20G; Murphy et al. 
2010) to extend these studies to a sample of nearby (z < 0.08) 
compact sources, which are selected at high frequency and so are 
expected to be the youngest and most recently-triggered radio 
AGN in the local Universe. This provides an important comple- 
ment to earlier studies, allowing us to improve our knowledge of 
the local population of associated H I absorption-line systems. 
Allison, Sadler, & Whiting (201 1) tested an automated spectral- 
line finding method, based on Bayesian inference and using 
simulated data from the Australian Square Kilometre Array 
Pathfinder (ASKAP; Deboer et al. 2009), as part of the prepa- 
ration for the First Large Absorption Survey in H I (FLASH) 1 . 
The data obtained here, from the Compact Array Broadband 
Backend (CABB; Wilson et al. 201 1) on the Australia Telescope 
Compact Array (ATCA), provide another test of this spectral- 
line finding method. We use these data to characterize any lim- 
itations that can arise from sources of systematic error, such as 
imperfect band-pass calibration and continuum subtraction. 

Throughout this paper we adopt a flat ACDM cosmology 
with fi M = 0.27, Q A = 0.73 and H = 71kms~ 1 Mpc" 1 . All 
uncertainties refer to the 68.3 per cent confidence interval, unless 
otherwise stated. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

2.1.1 Sample selection 

The sources observed in this programme were selected at 
20 GHz from the AT20G survey catalogue (Murphy et al. 2010), 
in order to target compact, core-dominated radio galaxies for 
which the 21 cm Hi line falls within the frequency range avail- 
able at the ATCA. The following selection criteria were used: 

(i) We used the list of AT20G sources identified by Sadler et 
al. (in preparation) with galaxies from the Third Data Release 
of 6dFGS Galaxy Survey (6dFGS DR3; Jones et al. 2009). This 
yielded a list of 202 galaxies south of declination 0° and with 
Galactic latitude |6| > 10°. 

(ii) We removed 42 sources from this list which were flagged 
in the AT20G source catalogue (Murphy et al. 2010) as having 
extended radio emission at 20 GHz (i.e. radio sources larger than 
10-15 arcsec in angular size). 

(iii) Finally, we restricted our target list to objects that lie in 
the redshift range 0.04 < z < 0.08, were at declination south of 
—20°, and had a low-frequency flux density of at least 50mJy 
in the National Radio Astronomy Observatory (NRAO) Very 
Large Array (VLA) Sky Survey (NVSS 1.4 GHz; Condon et al. 
1998) or the Sydney University Molonglo Sky Survey (SUMSS 
843 MHz; Mauch et al. 2003). This gave a final candidate list of 
45 objects. 

Of these 45 candidate sources, three have already been 
searched for H I absorption by Morganti et al. (2001) and one by 
van Gorkom et al. (1989). Eight of the 45 are extended sources in 
the lower-frequency NVSS or SUMSS images, and these were 
given lower priority, since our focus in this project is on com- 
pact sources. Table 1 gives a list of properties for these 45 can- 
didate sources, of which 29 were observed using the ATCA in 
2011. Note that two of the sources observed in our programme 
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have extended structures at low frequency (J 1 3 1 1 24^-42240 and 
J23 1905^120648 are classified as FR-I). 



2.1.2 Observations with the ATCA 

Observations of the 29 target sources were carried out in 2x 
three contiguous 24 hr periods over 2011 February 03 - 06 and 
2011 April 23 - 26, using the ATCA. The targets were cho- 
sen to give preference to objects with compact radio emission 
at both 20 GHz and 1 GHz, while for the sake of efficient ob- 
serving also retaining a roughly even distribution in right ascen- 
sion. We observed using the CABB system, which provided a 
64 MHz zoom band across 2049 channels centred at 1 .342 GHz, 
equating to a velocity resolution of ~ 7km s _1 . This set-up al- 
lowed us to simultaneously search for H I absorption over the 
redshift range 0.0338 < z < 0.0844 with a single observ- 
ing band. The 6-element Compact Array was arranged in the 
6A East- West configuration, with baseline distances in the range 
0.34 - 5.94 km. At v = 1.342 GHz, this corresponds to an 
angular-scale sensitivity range of approximately 8-35 arcsec, 
and the 22 m antennas provide a primary beam full- width at half- 
maximum (FWHM) of 35 arcmin. 

The 29 target sources were split into two sub-samples of 
15 objects and observed in February and April respectively. 
The previously known associated absorption system J 1 8 1 934— 
634548 (also known as PKS 1814-637) was observed in both 
sub-samples, as a reliability check of the system and data reduc- 
tion pipeline. The target sources were observed in lOmin scans 
and interleaved with 1 .5 min observations of a nearby bright cali- 
brator for gain and band-pass calibration. The primary calibrator, 
PKS 1934-638 (5i. 38 4GHz = 14.94 ± 0.01 Jy) 2 , was observed 
at regular intervals (approximately every 1.5 - 2 h) for calibra- 
tion of the absolute flux scale, and also as a secondary calibrator. 
Each target source observation had a total integration time of ap- 
proximately 3 - 4 h (depending on the duration of local sidereal 
time for which the object was in the sky) producing an average 
noise-level per spectral channel of approximately 4 mjy. 



2.2 Data reduction 

2.2.1 Flagging and calibration 

Flagging, calibration and imaging of the data were performed 
using tasks from the MlRIAD 3 package (Sault et al. 1995), and 
implemented using a purpose-built CABB Hi data reduction 
pipeline written in PYTHON. The CABB data were simultane- 
ously observed in two 2049 channel frequency bands (contin- 
uum 1.075 - 3.124 GHz and zoom 1.310 - 1.374 GHz), with sys- 
tem temperatures only recorded for the continuum band. In order 
to ensure reliable on-the-fly correction of the zoom band visibil- 
ities we apply the system temperature measurements recorded 
in the continuum band using the MlRIAD task ATTSYS. To first 
order this correction is assumed to be representative of the true 
variation in system temperature since the continuum band in- 
cludes the entire zoom band within its frequency range. 

In order to ensure good calibration of the antenna gains 
and band-pass, we manually flag all source scans that do not 
contain an associated scan of a calibration source. We also flag 
CABB channels that are known to contain significant radio- 
frequency interference (RFI). Transient RFI signals are identi- 
fied and flagged using the task MIRFLAG, where the thresh- 
old is set at a reasonable level so as to avoid accidental removal 



2 http://www.narrabri.atnf.csiro.au/calibrators/ 

3 www.atnf.csiro.au/computing/software/miriad7 



© 2012 RAS, MNRAS 000, 1-16 



Hi absorption in AT20G compact radio galaxies 3 



Table 1. Properties of candidate AT20G sources for associated Hi absorption, within the redshift range 0.04 < z < 0.08. S20 is the AT20G 20 GHz flux 
density (Murphy et al. 2010). S 1A is the NVSS 1.4GHz flux density (Condon et al. 1998). So. 843 is the SUMSS 843 MHz flux density (Mauch et al. 2003). 
The corresponding spectral indices between the two lower-frequencies and 20 GHz are given by af \ and aQ° 843 respectively. z opt is the optical spectroscopic 
redshift. SC is the spectral classification as defined by Mahony et al. (2011). is the Two Micron All Sky Survey (2MASS) ii"-band absolute magnitude 
(Skrutskie et al. 2006). The 20 GHz compactness (6km v i s ) is defined as the ratio of the measured visibility amplitude on the longest ATCA baselines to that 
on the shortest baselines (Massardi et al. 201 1). FR is the Fanaroff-Riley type at low frequency, as defined by Fanaroff & Riley (1974); Baldi & Capetti (2009); 
Ghisellini (201 1). 
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J132920-264022 


13 29 20.70 


-26 40 22.2 


101 


102.4 




-0.01 




0.0502° 


Aa 


-25.14 


1.03 





J15 1741-242220 *% 


15 17 41.76 


-24 22 20.2 


3449 


2041.9 




0.20 




0.0490° 


Aae 


-25.71 


1.00 





J164416-771548 *t 


16 44 16.03 


-77 15 48.5 


399 




1165.9 




-0.34 


0.0427° 


Aae 


-25.06 


1.00 


11 
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* Source searched for H I absorption by van Gorkom et al. (1989) 
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of astrophysical spectral lines. This automatic RFI flagging was 
performed both before and after calibration. 

We calculate time-varying band-pass and gain corrections 
for the amplitude and phase of each calibrator. In additon 
gain solutions were also calculated for the bright target source 
Jl 8 1934-634548 (Si. 38 4 GHz = 13.3 ± 0.3 Jy) 2 . The absolute 
flux scale and spectral index corrections were applied to each 
secondary calibrator, based on the flux model for PKS 1934- 



638. Finally the calibration solutions for each secondary cali- 
brator were applied to the paired target sources. 



2.2.2 Spectral ripple subtraction 

In order to measure the optical depth of 21 cm H I absorption we 
need to extract a spectrum at the source position that contains 
a reliable estimate of the continuum and spectral-line compo- 
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nents. However the ATCA has a relatively large field of view 
and wide spectral band and so we often find additional signal in 
the extracted spectrum. This is caused by the presence of other 
bright sources within the field, which are convolved with the 
synthesised beam of the telescope, and so contribute frequency 
dependent side-lobes at the position of the target source. These 
spectral-ripple features will generate confusion with any broad, 
low optical depth spectral lines, and so introduce systematic er- 
ror in estimates of the continuum flux. 

To remove the spectral ripple we construct a multi- 
frequency synthesis image, with width approximately 6 times 
the primary beam (FWHM), of the brightest continuum com- 
ponents by using the Mfs option in the tasks INVERT and 
CLEAN. The CLEAN algorithm (Hogbom 1974) was performed 
on the non-deconvolved continuum image within 60 arcsecond 
boxes (using a cutoff based on an estimate of the rms in the 
image) at the positions of known sources from the SUMSS 
(5843 MHz > 8 mJy; Mauch et al. 2003) and NVSS (5*1.4 GHz > 
2.5 mjy; Condon et al. 1998) catalogues. We then mask the re- 
sulting continuum model within a 30 arcsecond radius of the tar- 
get source centroid (approximately 6 times the FWHM of the 
synthesised beam), using the tasks IMMASK and IMBLR, and 
subtract this model from the visibility data using UVMODEL. As 
an example, Fig. 1 shows the visibility amplitude spectrum of 
J223931-360912, averaged over all epochs and telescope base- 
lines, before and after applying this subtraction procedure. Since 
this process does not produce a perfect subtraction of the ripple 
component we expect some residual signal to remain at approx- 
imately 1-2 per cent of the target continuum signal. This will be 
in addition to any error generated from imperfect calibration of 
the band-pass. 



2.2.3 Imaging 

We construct data cubes from the flagged, calibrated and ripple- 
subtracted visibility data. We apply natural weighting to the vis- 
ibilities in order to optimise the signal-to-noise in the image 
data. Each cube was constructed with 2 arcsec pixel resolution, 
in order to provide good sampling of the synthesised beam. The 
image data were deconvolved from the synthesised beam using 
the CLEAN algorithm, to produce images with a restoring beam 
FWHM of 10 arcsec, typical of the maximum baseline distance. 
The target sources are mostly unresolved and so we extract the 
synthesised beam-weighted spectra from the 1 1 central pixels of 
each image plane. Fig. 2 displays 800 channels (equivalent to a 
velocity range of approximately 5600 km s -1 ) of the extracted 
spectrum for each of the target sources, centred on the expected 
position of the 21 cm H I absorption line from optical redshift es- 
timates. For the subsequent quantitative analysis presented here 
the central 1649 channels of the CABB zoom band are used, 
where 200 channels are flagged at each end of the band, to ac- 
count for possible edge effects. 

We estimate the standard deviation per spectral channel us- 
ing the median of the absolute deviations from the median value 
(MADFM). This statistic is a more robust estimator of the true 
standard deviation than the rms when a strong signal is present 
in the data. Under the assumption of normal distributed data, the 
MADFM statistic (p) is related to the standard deviation (it) by 
(Whiting 2012) 



v (MHz) 

1360 1350 1340 1330 1320 



p = \/2erf _1 (0.5)cr w 0.674cr, 



(1) 



and hence we can estimate the standard deviation in the spectral 
data by multiplying the MADFM statistic by a factor of 1.483. 
We calculate this statistic over the spectrum where we have sub- 
tracted the continuum and removed those channels that are ex- 
pected to contain the spectral line. We assume that the standard 
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Figure 1. The visibility amplitude versus the Doppler corrected bary cen- 
tric velocity, averaged over all epochs and array baselines, from the 
ATCA observations of J223931-360921. Top: Before removal of the 
spectral ripple. Bottom: After removal of the spectral ripple by subtract- 
ing a continuum model of other bright sources within the field of view. 
The spectral frequency is indicated on the upper abscissa; however it 
should be noted that the scale is linear only in velocity. 



deviation per channel is constant for the CABB data. Wilson 
et al. (201 1) show that there is minimal overlap in the signal be- 
tween the CABB spectral channels and so we assume that the 
uncertainty per channel is essentially uncorrelated in the subse- 
quent analysis. 



3 ANALYSIS 

3.1 Bayesian inference 

In the following analysis we adapt the method presented by Al- 
lison et al. (201 1) in order to account for the relatively complex 
line-profiles expected in real spectral-line data from associated 
H I absorption. Analytical models are assumed for both the con- 
tinuum and spectral-line components and simultaneously fit to 
the data using Bayesian inference. 

The posterior (or joint) probability for a set of model pa- 
rameters (9), given the data (d) and the model hypothesis (M), 
can be calculated using Bayes' Theorem, 



Pi(e\d,M) = 



Pr(d\0,M)Pr(0\M) 
Pr{d\M) 



(2) 
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Figure 2. Synthesised beam-weighted spectra extracted from the position of each target source. For each spectrum the ordinate shows the synthesised 
beam-weighted flux (in Jy) and the abscissa the Doppler corrected barycentric velocity (in 10 4 km s — 1 ). For clarity, only 800 of the 2049 CABB zoom 
band channels are displayed, centred near the optical redshift given in Table 1 (vertical dashed line). The grey region represents an estimate of the 1 a 
uncertainty in the optical redshift. 
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The probability of the data given the set of model parameters, 
known as the likelihood, can be calculated based on assumptions 
about the distribution of the uncertainty in the data. If the data 
set is large and therefore quasi-continuous, one can approximate 
the likelihood by the form given for normal multivariate data 
(see e.g. Sivia 2006), 

L = Vx(d\0,M) 
1 



\/NT 



exp 



(d-myC- L (d-m) 



(3) 



where N is equal to the length of d, C is the covariance matrix 
of the data, C | is the determinant of the covariance matrix, and 
m is the vector of model data. In the special case where the un- 
certainty O) per datum is constant and uncorrelated, the above 
expression for the likelihood reduces to 



L = 



1 



: exp 



1 



£(4 



(4) 



The probability of the model parameters given the hypothesis, 
Pr(0\M), is known as the prior and encodes a priori informa- 
tion about the parameter values. For example, consider the situa- 
tion where the frequency position of an intervening H I absorber 
has been relatively well constrained from previous observations. 
If we trust these observations we might then choose to apply a 
normal prior to the spectral line position based on the known 
level of uncertainty. We would otherwise apply uninformative 
priors to the parameters if we were previously unaware of their 
value. Uninformative priors are typically uniform in either linear 
space (for location parameters) or logarithmic space (for scale 
parameters, known as the Jeffreys prior). 

The normalisation of the posterior probability in Equation 2 
is equal to the probability of the data given the model hypothe- 
sis and is referred to throughout this work as the evidence. The 
evidence is calculated by marginalising the product of the likeli- 
hood and prior distributions over the model parameters. This is 
given by 

E = Pr(d\M) 

= [ Pr(d\0,M)Pr(0\M)d9, 



(5) 



which follows from the relation given by Equation 2 and that 
the integrated posterior is normalised to unity. When the model 
hypothesis provides a good fit to the data the likelihood peak 
will have a large value, and hence the model hypothesis will 
have a large associated evidence value. However if the model is 
over-complex then there will be large regions of low likelihood 
within the prior volume, thus reducing the evidence value for 
this model, in agreement with Occam's razor. Estimation of the 
evidence is often key in providing a tool for selection between 
competing models. 

3.2 Application to spectral data 

We wish to investigate whether the spectral data warrant the 
inclusion of the spectral-line model in addition to the contin- 
uum model, and if so how many spectral-line components are 
favoured. The analytical expression for the continuum model is 
given by an n po i y 'th order polynomial of the form 



= S 



1 + 



n poly , V q 

?> 5-0 

I — 1 v ' 



(6) 



where S CO nt is the model continuum flux density as a function 
of the Doppler corrected barycentric velocity v, vo is the fidu- 
cial velocity given by the first spectral channel (equal to approx- 
imately 1.16 x 10 4 kms _1 for the data presented here), So is 



Table 2. Priors used in spectral-line detection. The model parameters are 
defined by Equations 6 and 7. (S bs) is the mean observed continuum 
flux density. 



Parameter 


Prior type 


Prior range 


So 


Uniform 


(0.1- 10) x (Sobs) 


si 


Uniform 


±1 


Vi 


Uniform 


min(u) - max(v) 


Au PWH M,i 


Uniform 


0.1- 10 4 kms- 1 


AS; 


Uniform 


(0 - 1) x So 


Systematic Error 


Normal 


1 ± 0.01 



the spectral flux density at vo and, Si is the coefficient for the ith 
term. For the purposes of the work presented here we choose to 
use a 1st order polynomial (i.e. linear) model for the continuum 
component of the spectrum. 

The spectral-line model is given by the sum over re comp 
Gaussian components of the form 



ASlinc 



A Si exp 

i=l 



-41n(2) 



(v 



(Al>FWHM,' 



(7) 



where AS; is the peak depth in flux density, «; is the velocity, 
and At>FWHM,i is the FWHM of the ith component. The set of 
model data (m) are calculated by taking the mean of the com- 
bined spectral model over each spectral channel. This assumes 
that the channel response of the CABB Filterbank Correlator is 
flat, which is a reasonable approximation for our purposes here 
(Wilson etal. 2011). 

We use uninformative priors for all of the parameters in 
our spectral model (see Table 2). The maximum value for the 
continuum normalisation So is set by an estimate of the mean 
flux of the source. The line-depth prior on ASi is set by the 
physical limit of the mean flux density of each source. The ve- 
locity position (Vi) is limited to within the range of velocities 
recorded in the data. The maximum velocity width (Aufwhih,;) 
of 10 4 fans -1 is considered to be a physically reasonable limit. 

We use the following relationship for the relative probabil- 
ity of two hypotheses, 



(8) 



Pr(.Mi|ri) _ Pr(d\Mi) Pr(Mi) _ gi Pr(Mi) 
PT(M 2 \d) ~ Pr(d|X 2 ) Pr(M 2 ) ~ E 2 Pr(_M 2 ) ' 

to quantify the relative significance of the combined spectral line 
& continuum model versus the continuum-only model. The ratio 
Pr(M i ) /Pr(A^2) encodes our prior belief that one hypothesis 
is favoured over another. Since we assume no prior information 
on the presence of spectral lines, this ratio is equal to unity, and 
so the above selection criterion is then equal to the ratio of the 
evidences. We define the quantity 



R= In 



Em 

E CQn t 



(9) 



where Em is the evidence for the combined spectral line & con- 
tinuum hypothesis and E cont is the evidence for the continuum- 
only hypothesis. A value of R greater than zero gives the level 
of significance for detection of the spectral line. A value of R 
less than zero indicates that the data do not warrant the inclusion 
of a spectral-line component and so the detection is rejected. 



3.3 Implementation 

We fit our spectral model using the existing MultiNest 4 pack- 
age developed by Feroz & Hobson (2008) and Feroz, Hobson, 



http://ccpforge.cse.rl.ac.uk/gf/project/multinest/ 
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Table 3. The inferred values for model parameters from fitting a multiple Gaussian spectral-line model, 
where n comp is the best-fitting number of components; Vi is the component velocity; Atj FW hm,i is the 
FWHM; ASi is the depth and _R is the natural logarithm of the ratio of probability for this model versus the 
no spectral-line model. Only target sources for which we detect H I absorption are shown. 



AT20G name 


^comp 


(kms -1 ) 


(kms -1 ) 


ASi 
CmJv) 


R 


JUO 1 605—234352 


2 


19197.4 ± l.o 


64.5 ± 4.4 


41.0 ± 2.8 


coco I n 1 

535.3 ±0.1 






19286 ± 14 


138 ± 23 


11.5 ± 1.5 




J181934-634548 a 


4 


19161.6 ±0.6 


46.0 ±0.5 


1812 ±20 


103157.2 ±0.2 






19173.8 ±2.4 


212.5 ± 8.3 


196.5 ±8.3 








19182.6 ±0.1 


19.4 ±0.4 


1335 ± 32 








19331.6 ± 2.2 


25.9 ± 7.5 


102 ± 12 




J181934-634548 6 


4 


19162.9 ±0.2 


47.1 ± 0.3 


1869 ± 15 


211968.3 ±0.2 






19179.1 ± 1.7 


224.1 ± 5.0 


175.0 ±4.9 








19183.4 ±0.1 


19.4 ±0.3 


1288 ± 22 








19333.3 ±0.8 


18.8 ±2.6 


107.8 ±9.2 




J20540 l^t-24238 


1 


13048 ± 34 


520 ± 140 


4.5 ±0.7 


26.3 ±0.1 



a Source observed in 201 1 February. 
b Source observed in 201 1 April. 



& Bridges (2009). This software uses nested sampling (Skilling 
2004) to explore parameter space and robustly calculate both the 
posterior probability distribution and the evidence, for a given 
likelihood function and prior (provided by the user). 

We first test for detection of the absorption line using a 
single component spectral-line model with multi-modal nested 
sampling, whereby samples are taken of multiple likelihood 
peaks within parameter space, thus allowing for the possibility 
of multiple absorption lines in the spectrum. For each peak in 
likelihood we calculate a local Em value for the single compo- 
nent model, and then compare it with _E C ont for the continuum- 
only model. The significance of the combined spectral line & 
continuum model is inferred by the relative local value of Em 
compared with -Econt- If the local value of Em is less than or 
equal to E cont then this detection is rejected. Following the suc- 
cessful completion of the nested sampling algorithm, both the 
multi-modal local evidence values and the model parameter pos- 
terior probability are recorded. 

For each detection obtained using the single component 
spectral-line model, we repeat the fitting process by iteratively 
increasing the number of components and comparing the values 
of R. This process continues until we optimize R and we have 
found the best-fitting number of Gaussian components for the 
spectral line. 

3.4 Systematic errors 

We estimate an upper limit of 1 per cent for the calibration error 
in the CABB spectral data (Wilson et al. 201 1). This is consistent 
with the difference in flux density (approximately 0.3 per cent) 
for February and April observations of Jl 8 1934-634548. We ac- 
count for systematic errors introduced by imperfect band-pass, 
gain and flux-scale correction by introducing an additional pa- 
rameter that multiplies the model spectral flux density. This nui- 
sance parameter has a normal prior distribution with a mean 
value of unity and a standard deviation equal to 1 percent, and is 
subsequently marginalised when calculating the posterior prob- 
ability distributions for the model parameters. 

3.5 Continuum residuals 

When searching for H I absorption in real interferometric data 
we expect to obtain some false detections due to residual signal 
in the continuum that is not well fit by our continuum model. 



This could be due to imperfect spectral-ripple subtraction (see 
Section 2.2.2) or band-pass calibration. We expect that this prob- 
lem will be worse for models with significantly large spectral- 
line widths and low depths, since at these values we would es- 
sentially be improving the continuum fit. The Bayesian method 
generates a marginalised posterior distribution for each model 
parameter, including the spectral-line width. We therefore in- 
spect the line- width parameter a posteriori for unphysical values 
and potential association with a broad continuum residual. Since 
we are looking for associated absorption we can also reject those 
candidate broad detections with a line velocity that is very dif- 
ferent to that predicted by existing optical redshift estimates. 

Lastly we invert the spectrum and search for emission-like 
features, which could arise either from emission lines or con- 
tinuum artefacts. We expect that a broad absorption-like ripple 
feature would likely be accompanied by an emission-like rip- 
ple feature. The residual contamination of spectra, both from 
spectral-ripple subtraction and imperfect band-pass calibration, 
will be of concern for future large field of view, wide-band H I 
absorption surveys 5 . These surveys will require a reliable sky 
model for continuum subtraction. 



4 RESULTS AND DISCUSSION 
4.1 Observational results 

Fig. 3 summarizes both the absorption and emission features that 
were detected in data from all 29 target sources, using a sin- 
gle Gaussian spectral-line model and multi-modal nested sam- 
pling. It is evident from this plot that there is a significant num- 
ber of identified false detections, clustered at broad spectral-line 
widths (Aijfwhm > 10 3 kms _1 ) and low peak optical depths 
(Tpoak < 0.1), as was discussed in Section 3.5. The probability 
distributions for some of these detections are peaked beyond the 
maximum prior for the spectral-line width (10 4 km s _1 ) and so 
these points are shown at the limiting value. We obtain a total of 
59 false detections, of which 30 are absorption-like and 29 are 
emission-like spectral features. 

The false detections shown in Fig. 3 form a distinct group- 
ing of broad, shallow features, which are clearly separate from 
the significant Hi absorption detections in J001605-234352 



5 http://www.physics.usyd.edu.au/sifa/Main/FLASH 
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Figure 3. Best-fitting (maximum likelihood) values for the peak ob- 
served optical depth versus spectral-line width for absorption-like (cir- 
cles) and emission-like (crosses) features detected using a single Gaus- 
sian model and multi-modal nested sampling. The filled circles are iden- 
tified as H I absorption, while empty circles represent false detections. 
The red points represent features detected from observations of the 
brightest continuum source J181934-634548. The dashed lines indicate 
the position of the single channel width of 7 km s — 1 and an optical depth 
of t = 1 respectively. 



and J181934-634548. The weaker detection seen in J205401- 
424238 has best-fitting values of peak observed optical depth 
and spectral-line width that are closer to those of the grouping of 
false detections. However it is the narrowest of the absorption- 
like features in this group (Aufwhm = 520 ± 140 km s^ 1 ) 
and has a velocity that is consistent with that expected from the 
optical spectroscopic redshift. In total we have detected three 
associated H I absorption lines, of which two were previously 
unknown, giving a 10 per cent detection rate. Fig. 4 shows the 
optical and radio continuum images, and spectral model fitting 
results for the three sources with detected associated H I absorp- 
tion. Table 3 summarizes the properties of the individual Gaus- 
sian components from model fitting. 

The optical depth is related to the continuum flux density 
(Scant) and spectral-line depth (ASimc) by 



t(v) 



■In 1 



fS coni 



(10) 



where / is the covering factor of the background source. 
Under the assumption that the absorption is optically thin 
(AS'iine/JScont < 0.3), the above expression for the optical 
depth reduces to 



t(v) 



fS con 



Tobs(f) 



(ID 



where t ^ b (v) is the observed optical depth as a function of ve- 
locity. The Hi column density (in cm -2 ) is related to the inte- 
grated optical depth (in kms -1 ) by (Wolfe & Burbidge 1975) 



JVhi = 1.823 x 10 18 T spin / r(v)dv 

,18 Tspin 



1.823 x 10 1 



T ob a(v)dv, 



(12) 



where T sp i n is the mean harmonic spin temperature of the gas 
(in K). Table 4 summarizes these quantities derived from model 
fitting. 



4.2 Notes on individual sources 

J001605-234352: We detect a previously unknown H I absorp- 
tion feature with a peak observed optical depth of 0.169 ± 
0.007 and an integrated optical depth of 16.9 ± 0.7 km s -1 . 
van Gorkom et al. (1989) include this source in their sample 
of radio galaxies with compact cores, observed using the Very 
Large Array. However, they do not include it in their search for 
21 cm Hi absorption since the unresolved continuum flux den- 
sity (255 mJy) did not meet their selection criterion of greater 
than 350 mJy. The absorption line feature is best fit by two 
Gaussians, with a narrow component (Aufwhm = 64.5 ± 
4.4 km s -1 ), located at the spectroscopic optical redshift, and a 
single broad component (A«fwhm = 138 ±23 km s _1 ), which 
is slightly redshifted by 96 ± 46 km s ~~ 1 . The host galaxy is clas- 
sified as a blue spiral (UK Schmidt Telescope, Automated Pho- 
tographic Measuring Bright Galaxy Catalogue; Loveday 1996) 
located within the galaxy cluster Abell 14 (Savage & Wall 1976; 
Dressier 1980). Fig. 4 shows the blue optical image at the posi- 
tion of the target source (UK Schmidt Telescope, SuperCosmos 
Survey; Hambly et al. 2001), with an edge-on spiral galaxy that 
exhibits some disruption at the edges. 

Both the H I absorption line profile and the host morphol- 
ogy for this source demonstrate strong similarity with J181934- 
634548, studied at multiple wavelengths by Morganti et al. 
(2011). The narrow, deep absorption component is consistent 
with a gas-rich galactic disc, while the redshifted broad shal- 
low component could be attributed to higher velocity gas close 
to the central AGN, perhaps in a circumnuclear disc. However to 
verify this interpretation we require further 21 cm observations 
at higher angular resolution (via long baseline interferometry). 

J003908-222001: We obtain a low significance spectral-line 
feature (R = 4.26 ± 0.11) at z = 0.0785 ± 0.0031 with a 
width of Aufwhm = 2.7 ±1.2 km s _1 . In Fig. 3 this feature is 
isolated to the very narrow spectral-line width and high optical 
depth region. This feature is redshifted by 4230 ± 930 km s -1 
with respect to the systemic velocity, given by the spectroscopic 
optical redshift (z = 0.06438 ± 0.00025; Vettolani et al. 1989), 
and, if attributed to H I absorption, would likely represent an in- 
fall of gas at very high velocity. Because of this, and that the 
spectral-line width spans a single spectral channel, we suspect 
that this feature arises due to relatively high noise in this chan- 
nel. 

J181934-634548: This compact steep spectrum source 
(Tzioumis et al. 2002) has a previously known strong absorption 
feature (Veron-Cetty et al. 2000; Morganti et al. 2001) and was 
studied extensively at multiple wavelengths and angular scales 
by Morganti et al. (201 1). In particular the host galaxy is known 
to have a strong edge-on disc morphology with heavily warped 
edges, consistent with recent merger activity. This is not easily 
apparent upon inspection of the blue optical image shown in 
Fig. 4, due to the presence of a bright coincident star. 

We re-detect the absoiption feature in both of our February 
and April observations, with derived peak and integrated optical 
depths that are consistent within the error. We estimate a peak 
observed optical depth of 0.193 ± 0.001 and an integrated opti- 
cal depth of 11.92 ± 0.05 km s . The absorption line feature is 
best fit by four Gaussians, consisting of a deep, narrow compo- 
nent (combination of two Gaussians; Aufwhm ~ 20kms~ 1 
and Aufwhm ~ 50kms _1 ), a shallow, broad component 
(Aufwhm ~ 200kms _1 ), and a small redshifted component 
(Avfwhm ~ 20kms _1 ). Of the two strong features, the broad 
component is slightly redshifted with respect to the narrow com- 
ponent, generating an asymmetry in the wings of the line profile 
towards the red end of the spectrum. Morganti et al. (201 1) used 
observations with the Australian Long Baseline Array (LB A) 
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Figure 4. Target sources with detected H I absorption. Left: Images from the ATCA observations (contours) and the SuperCosmos Sky Survey (grey- 
scale; Hambly et al. 2001). The contours represent intervals of lOpercent of the peak continuum brightness. The FWHM of the ATCA synthesised 
beam is displayed in the bottom left-hand corner of each image. Right: The corresponding spectra from the ATCA observations. The data have been 
simultaneously fit by a combined continuum and multiple Gaussian spectral-line model. The solid black line represents the data after subtraction of the 
best-fitting continuum component. The solid red line represents the best-fitting multiple Gaussian spectral-line model, and the dashed red lines represent 
individual components. The solid blue line represents the best-fitting residual, including a velocity-axis offset for clarity. R refers to the detection 
significance as defined by Equation 9. 
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Table 4. A summary of derived H I absorption properties from model fitting, where (J chan is the estimated uncertainty 
per channel; S cont> p Ca k is the spectral flux density of the continuum model at either the position of peak absorption 
or the optical redshift estimate; AS pcak is the peak spectral-line depth; r obs pcak is the observed peak optical depth 
and J r b s di) is the observed velocity integrated optical depth. Upper limits are calculated assuming a single Gaussian 
spectral line of FWHM equal to 30 km s ~~ 1 and peak depth equal to 3 cr chan . 
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J131 124-^42240 


3.28 


75.9 ±0.8 


<9.8 


< 0.130 


<4.13 


< 18.88 


J132920-264022 


3.23 


73.5 ± 0.7 


<9.7 


< 0.132 


<4.21 


< 18.88 


J171522-652018 


3.26 


150.7 ± 1.5 


<9.8 


< 0.065 


< 2.07 


< 18.58 


J180957^455241 


3.40 


786.9 ± 7.0 


< 10.2 


< 0.013 


< 0.41 


< 17.88 


J181934-634548 a 


12.55 c 


13709 ± 57 


2632 ± 15 


0.192 ±0.001 


11.92 ±0.06 


19.337 ± 0.002 


J181934-634548 & 


8.87 c 


13754 ±55 


2654 ± 13 


0.193 ±0.001 


11.93 ±0.04 


19.337 ± 0.001 


J192043-383107 


3.37 


199.0 ± 1.9 


< 10.1 


< 0.051 


< 1.62 


< 18.47 


J205401^24238 


3.47 


162.1 ± 1.5 


4.5 ± 0.7 


0.028 ± 0.004 


13.4 ±2.0 


19.38 ± 0.07 


J205 837-575636 


3.40 


221.2 ±2.1 


< 10.2 


< 0.046 


< 1.47 


< 18.43 


J210602-712218 


3.38 


1021.9 ±9.0 


< 10.2 


< 0.010 


< 0.32 


< 17.76 


J220253-563543 


3.23 


51.8 ±0.5 


<9.7 


< 0.187 


<5.98 


< 19.04 


J223931-360912 


4.22 


377.8 ±3.6 


< 12.7 


< 0.034 


< 1.07 


< 18.29 


J23 1905-420648 


3.42 


270.7 ± 2.6 


< 10.3 


< 0.038 


< 1.21 


< 18.34 


J233355-234340 


4.47 


606.2 ± 5.7 


< 13.4 


< 0.024 


< 0.71 


< 18.11 


J234129-291915 


4.11 


204.9 ± 2.0 


< 12.3 


< 0.060 


< 1.92 


< 18.54 



a Source observed in 201 1 February. 
b Source observed in 2011 April. 

c Significantly higher per-channel variance due to band-pass calibrator. 



to demonstrate that the absorption consists of a deep, narrow 
component from the edge-on disc of the galaxy, and a spatially 
resolved shallow, broad component from the cirumnuclear disc 
close to the radio source. 

In addition to these principle features we also detect a 
small component, which is redshifted from the systemic veloc- 
ity by approximately 100 km s -1 , and has a peak optical depth 
w 0.008 and integrated optical depth w 0.2kms _1 . The data 
from both of our February and April observations favour the in- 
clusion of this additional Gaussian component, with the natural 
logarithm of the ratio of the probabilities of a 4-component to a 
3-component model of 36.1±0.2 and 83.1±0.2 respectively. In- 
spection of the LBA spectrum obtained by Morganti et al. (201 1) 
shows a similar feature, at this velocity, in the broad red-shifted 
wing of the profile seen against the northern lobe. It is therefore 
likely that this feature represents non-uniformity in the gas on 
the red-shifted side of the circumnuclear disc. A less likely in- 
terpretation, since the feature is only seen against the northern 
lobe, is that it indicates the presence of infall of gas towards the 
central AGN at a velocity greater than 100 km s _1 . 

J205401-424238: We detect a previously unknown H I absorp- 
tion feature with a peak optical depth of 0.028 ± 0.004 and 
an integrated optical depth of 13.4 ± 2.0 km s -1 . This is the 
least significant of our three detections, however the single- 
component Gaussian fitting returns a significant relative log- 



probability of 26.3 ± 0.1 for the spectral line hypothesis, and 
has a position that is consistent with the spectroscopic opti- 
cal redshift. We note that this feature is apparent throughout 
the observation, in both polarisation feeds, and does not ap- 
pear in spectra that are extracted from spatial positions signifi- 
cantly off the source. The single component is shallow and broad 
(Awpwhm = 520 ± 140 km s -1 ) producing a relatively low 
peak, but high integrated, optical depth. 

The radio source is signified by a bright, flat spectrum, and 
was observed in the Combined Radio All-Sky Targeted Eight 
GHz Survey (CRATES; Healey et al. 2007). The host galaxy is 
classified as an elliptical by Loveday (1996), which is consistent 
with the optical image from the UK Schmidt Telescope shown 
in Fig. 4. That we detect only a single broad, shallow component 
at 21 cm suggests that we might be seeing only the absorption 
from gas close to the AGN, and nothing from the gas-poor ellip- 
tical host. Again in order to verify this interpretation we require 
further 21 cm observations at higher angular resolution (via long 
baseline interferometry). 



4.3 Stacking for a statistical detection 

While we do not detect H I absorption in 26 of the target sources 
with the ATCA, we can search for a statistical detection by stack- 
ing the individual spectra. For each spectrum we calculate the 
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Figure 5. The binned optical depth spectra for the 26 sources with non-detections of H I absorption. The binned velocity axis is centred on the rest 
velocity given by the optical redshifts (see Table 1). A velocity range of ±1500 kms - 1 is shown for clarity. Left: The data have been binned to 
10kms _1 , and so similar to the spectral resolution of the individual spectra. Right: The data have been binned to a resolution of 122 km s -1 , equal to 
the mean error in the optical spectroscopic redshift. The errorbars represent 1 a. 



observed optical depth spectra, using the the best-fitting con- 
tinuum model, and centre the velocity axis on the rest velocity, 
given by the optical spectroscopic redshift. We define the ob- 
served optical depth and weight, for the j th datum, by 



= 1 - 



^cont^' 
((Tj/ Scont,j) 



(13) 
(14) 



where Sj is the spectral flux density, S CO nt,j is the best-fitting 
model continuum flux density and Oj is the uncertainty in Sj . 
The binned optical depth (r hin ), velocity («bi n ) and uncertainty 
(obin) are then given by 



Tbin = 



^bin = 



Cbii 



1/2 



.3=1 



£' 

■3 = 1 



£' 

U=i 



(15) 
(16) 
(17) 



where Vj and z op t,j are the spectral velocity and optical spec- 
troscopic redshift respectively. Fig. 5 shows the binned flux den- 
sity versus velocity for the stacked non-detections. We consider 
two values for the velocity bin widths, in one case similar to 
the spectral resolution of an individual spectrum and the second 
case equal to the mean uncertainty in the optical spectroscopic 
redshift (« 122 km s -1 ), and hence approximately equal to the 
smearing scale for any spectral line. The optical depth sensi- 
tivity for the narrow-binned spectrum is 0.14 per cent. Using the 
spectral-line finding method outlined in Section 3, with a normal 
prior on the spectral-line velocity of v = ± 122 kms -1 , we 
find that there is no evidence of a statistical detection of 21 cm 
associated H I absorption in the stacked non-detections. 



4.4 Intervening H I absorption 

Based on both the lower redshift limit available from our ob- 
servations (z m in = 0.0338) and the optical spectroscopic red- 
shifts of the 29 target sources, the total redshift path searched 
for potential intervening H I absorption is Az to tai ~ 0.6093, 
assuming that there is no high-velocity, redshifted intervening 



gas towards our sources. We can estimate the expected number 
of high column density intervening H I absorbers in this redshift 
path from previous work on nearby galaxies. Zwaan et al. (2005) 
used 21 cm line observations of 355 nearby galaxies in the West- 
erbork H I survey of Spiral and Irregular galaxies (WHISP; van 
der Hulst, van Albada, & Sancisi 2001) to estimate the num- 
ber density of damped Lyman-a galaxies (DLAs, defined by 
Nm > 2 x 10 20 cm~ 2 ) at z = 0. Based on their estimate 
of dN/dz — 0.045 ± 0.006, and assuming no strong evolu- 
tion over our redshift range, we would expect to find approxi- 
mately 0.027 DLAs within our total redshift search path (assum- 
ing no a priori information of the location of foreground galax- 
ies). If we assume that T sp i n // = 100 K, then a single DLA 
equates to an observed integrated optical depth of 1.09kms -1 , 
which is comparable to the mean sensitivity of our sample of 
1.85 kms -1 . It is therefore highly unlikely that we would have 
detected intervening H I absorption in our sample, and that we 
could have confused one of these systems with a continuum 
residual signal (which would also require a significantly broad 
spectral-line width for the absorption). It should also be noted 
that this result does not change significantly if we use a dif- 
ferent estimate for the low redshift number density, for exam- 
ple dN/dz = 0.046±o.o2 (Ryan-Weber, Webster, & Staveley- 
Smith 2003) or dN/dz = 0.026 ± 0.003 (Braun 2012). 



4.5 Detection rate for associated H I absorption 

4.5.1 Comparison with other surveys 

Table 5 compares our results with those of other similar asso- 
ciated H I absorption line studies that preferentially target com- 
pact objects. The total detection rate for our survey (lOper cent) 
is the lowest out of these studies, and comparable to the survey 
of compact cores of E and SO galaxies by van Gorkom et al. 
(1989). Fig. 6 shows the observed integrated optical depth ver- 
sus the 1 .4 GHz flux density, for all the associated H I absorption 
surveys listed in Table 5 and the z > 0.1 sample discussed by 
Curran et al. (201 1). From this plot it is evident that the majority 
of our sample have relatively low fluxes and so our observations 
are typically sensitive to a higher observed integrated optical 
depth limit than the other surveys. This is a possible reason for 
the relatively low detection rate obtained in this work, indeed if 
one only counts detections with an observed peak optical depth 
greater than 0.05 then this number is comparable to the other sur- 
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Table 5. A comparison with other similar surveys of associated H I absorption. N \> s is the total number of target sources; N^ et is the number of H I absorption 
detections and iVo.05 is the number of detections with observed peak optical depth greater than 0.05. 



Ref. Telescope z range JV b s Ndet Rale Nq.05 Targets 

(per cent) 



van Gorkom et al. 1989 


VLA 


0.004 


-0.128 


29 


4 


14 


1 


Compact cores of E/S0 galaxies 


Morganti et al. 2001 


VLA/ATCA 


0.011 


- 0.220 


23 


5 


22 


2 


Radio galaxies selected at 2.7 GHz 


Vermeulen et al. 2003 


WSRT 


0.198 


- 0.842 


57 


19 


33 


2 


Compact radio sources 


Orienti et al. 2006 


WSRT 


0.021 


-0.668 


6 


2 


33 


1 


High-frequency peaker radio galaxies 


Gupta et al. 2006 


Arecibo/GMRT 


0.016 


-1.311 


27 


6 


22 


1 


GPS and CSS radio sources 


Emonts et al. 2010 


VLAAVSRT 


0.002 


- 0.044 


23 


6 


26 


3 


Low luminosity radio sources (compact and FR-I) 


Chandola et al. 2011 


GMRT 


0.024 


-0.152 


18 


7 


39 


3 


GPS and CSS radio sources (CORALZ) 


This work 


ATCA 


0.041 


- 0.076 


29 


3 


10 


2 


AT20G radio sources selected at 20 GHz 



veys (albeit with large uncertainty due to low-number statistics). 
However there are potential errors associated with this interpre- 
tation, since some of our observations are sensitive to observed 
optical depths much less than 0.05 (see Table 4) and there are 
some high integrated optical depth detections seen in the lower 
flux sources (Fig. 6). We also note that the observed optical depth 
upper limits are based on simple assumptions about the underly- 
ing spectral-line width and depth. 



4.5.2 Ultra-violet luminosity 

A low detection rate of associated H I absorption in high red- 
shift sources was found by Curran et al. (2008), who obtained 
no detections from the 13 (z > 3) radio galaxies and quasars 
searched. Upon a detailed analysis of the photometry of the 
objects, all of the targets were found to have ultra-violet lu- 
minosities greater than an apparently critical value (L1216 ~ 
10 23 WHz -1 ) at A = 1216 A, above which the photon energy 
is sufficient to excite the hydrogen atom above the ground state 
and close to ionisation 6 . Curran et al. (2008) therefore attributed 
the low detection rate to the high redshift selection of the targets 
biasing towards the most ultra-violet luminous objects. 7 In fur- 
ther work by Curran et al. (201 1) it was found that there are no 
detections of Hi absorption in the known 19 target sources with 
ii2is > 10 23 WHz -1 . Assuming a 50 per cent detection rate 
(see Curran & Whiting 2010) this distribution has a binomial 
probability of 9.5 x 10 -7 of occuring by chance. 

Fig. 7 shows the estimated ultra-violet luminosity, at a rest 
wavelength of A = 1216 A, for all the associated H I absoiption 
surveys. Most of our target sample (24 sources) have measured 
fluxes from the Galaxy Evolution Explorer (GALEX; Martin 
et al. 2003), with the remaining each having several flux mea- 
surements at rest-frame frequencies of v > 10 14 Hz, thus allow- 
ing us to determine the luminosity at 2.47 x 10 15 Hz. However 
only one of the sources has L1216 > 10 22 WHz" 1 ; J171522- 
652018 with log 10 (Li2i6) ~ 22.7. This is one of the three 
sources for which there is no photometry at v > 10 15 Hz and 
so the estimated ultra-violet luminosity is unreliable. Therefore 
the low detection rate of our survey cannot be explained by the 
same extreme ultra-violet luminosities introduced by selecting a 
high redshift sample. 



6 Curran & Whiting (2012) show that this is the critical luminosity re- 
quired to ionise all of the gas in a large galaxy. 

' Despite shortlisting the most optically faint objects (those with blue 
magnitudes of B > 19, see fig. 5 of Curran, Whiting, & Webb 2009), in 
the Parkes Quarter- Jansky Flat-spectrum Sample (PQFS; Jackson et al. 
2002). 




0.01 0.1 1 10 

Rest— frame 1 .4 GHz flux density [Jy] 

Figure 6. The normalised H I column density (1.823 X 10 18 / T b s d?;) 
versus the rest-frame 1.4 GHz flux density for the radio galaxies and 
quasars searched in 21cm associated Hi absorption. The filled sym- 
bols represent the 21 cm detections and the unfilled symbols the non- 
detections. The squares show results from previous searches (Table 5 
and Curran et al. 201 1) with the circles designating the sources presented 
here. 

4.5.3 Radio luminosity 

Fig. 7 also shows the estimated radio luminosity, at a rest wave- 
length of A = 21 cm, for all the associated Hi absorption sur- 
veys. Not suprisingly, given target source selection, there is an 
apparent strong correlation between the radio luminosity and the 
redshift (17 a). It is evident from this plot that the total detec- 
tion rate for H I absorption is higher at L21 cm > 10 26 WHz -1 
(30.8percent, versus 19.2percent at lower luminosities). The 
target sources in our sample have radio luminosities of L21 cm < 
10 26 WHz -1 , and so the lower detection rate in our sample 
might be consistent with that expected for these sources. We 
note that the data shown in Fig. 7 also include observations by 
Emonts et al. (2010), who conducted a survey of nearby low lu- 
minosity radio galaxies (22 < log 10 (Z/2icm) < 25WHz -1 ), 
and obtained a relatively high detection rate of 6 (26 per cent) 
out of a sample of 23 sources 8 . However, because all of these 
Hi absorption systems were found to be associated with the 
host galaxy's large-scale interstellar medium (ISM), their sam- 
ple may be biased against intrinsically weak sources that are 



We have not included the three tentative detections from this survey. 
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Figure 7. The rest-frame A = 1216 A continuum luminosity (top) and 21cm continuum luminosity (bottom) versus the host redshift for the radio 
galaxies and quasars searched in 21 cm associated H I absorption. The filled symbols/hatched histogram represent the 21 cm detections and the unfilled 
symbols/unfilled histogram the non-detections. The squares show results from previous searches (Table 5 and Curran et al. 2011) with the circles 
designating the sources presented here. 



confined by the surrounding gas. Given the selection effects and 
small number of targets associated with existing HI absorption 
surveys, it is difficult to compare and interpret their individual 
detection rates. For this, future proposed large-scale surveys, 
such as FLASH (see Section 1) will be crucial. 

Given the apparent higher total detection rate for H I ab- 
sorption in higher radio luminosity sources, we might speculate 
that there is a correlation between the neutral hydrogen column 
density of the host and the radio luminosity of the source. For ex- 
ample Chandola et al. (201 1) searched for a correlation between 
the neutral hydrogen column density and 5 GHz radio luminos- 
ity, both in the Compact Radio sources at Low Redshift sample 
(CORALZ; Snellen et al. 2004) and the survey by Gupta et al. 
(2006), and found no statistically significant evidence. However, 
since the neutral hydrogen column density is derived from an as- 
sumed T a pi n //, it is possible that higher spin temperatures in the 
most radio luminous sources nullify any underlying relationship. 
If so then we might expect to see a limiting radio luminosity 
above which no detections are obtained, similar to that seen for 
the ultra-violet luminosity. This was tested by Curran & Whit- 
ing (2010) who found no significant evidence for such a limiting 
radio luminosity value. 



4.5.4 Compact sources 

When selecting our AT20G candidate sources for observation, 
priority was given to objects that did not exhibit extended 
emission, both at 20 GHz and in the lower frequency NVSS 
(1.4 GHz) and SUMSS (843 MHz) images. The motivation to 
preferentially select young and compact sources was based on 
previous work, which claims a significantly higher detection 
rate in compact sources (see Section 1). Given these results we 
might conclude that the compact sources either have intrinsi- 



cally higher neutral hydrogen column densities (e.g. Pihlstrom 
et al. 2003) or that the flux density of the background contin- 
uum source is concentrated behind the absorbing gas, boosting 
the covering factor /. 



However Curran & Whiting (2010) showed that when the 
> 10 23 WHz _1 sources were removed from the sam- 
ple, the 21 cm detection rate in compact sources was not signif- 
icantly higher than that for non-compact sources. The elevated 
detection rate in compact sources was therefore interpreted as 
these having apparently lower ultra-violet luminosities, consis- 
tent with the hypothesis that these are young objects. In Fig. 8 
we update the results of Curran & Whiting (2010) to include the 
work presented here and more recent results (including Emonts 
et al. 2010; Curran et al. 2011; Chandola et al. 2011). This fig- 
ure shows the total number of detections and non-detections for 
sources classified as compact (either compact symmetric ob- 
ject - CSO, compact steep spectrum - CSS, gigahertz peaked 
spectrum - GPS or high frequency peaker - HPF), compared 
with those that are non-compact or unclassified. While the to- 
tal number of detections is greatest for the compact sources we 
note that there is also a significant number of non-detections. 
The detection rates for each classification, for the whole sam- 
ple, are 32.7percent (compact) versus 18.6percent (other) and 
16.1 percent (unclassified). For the sub-sample of sources with 
ii2i6 ^ 10 23 WHz _1 , the detection rates are 33.9 per cent 
(compact) versus 23.5 per cent (other) and 18.2 per cent (unclas- 
sified). Therefore the detection rate is slighty higher for the com- 
pact sources, although not significantly so if we only consider 
the sub-sample of sources with L1216 ^ 10 23 WHz -1 . 
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Figure 8. The incidence of detections (hatched histogram) and non- 
detections (unfilled histogram) for compact objects (CSO, CSS, GPS 
and HFP), compared to the others (non-compact and unclassified). Top: 
the whole sample (Table 5 and Curran et al. 2011). Bottom: those with 
ii2i6 «: 10 23 WHz- 1 . 



In particular it seems that for non-ideal band-pass calibration 
and continuum subtraction there will exist in the spectral data 
continuum residuals that look like broad, shallow spectral lines. 
For targeted associated absorption studies we are fortunate that 
we have a priori information on the expected spectral-line veloc- 
ity from measured optical redshifts. However in the case of fu- 
ture large-scale blind surveys for both intervening and associated 
absorption it will be vitally important to have good calibration 
of the band-pass and simultaneous accurate models of the con- 
tinuum sky. We note that by studying the posterior distribution 
of the model parameters, produced by the Bayesian inference 
method, we can filter out those detections that do not meet pre- 
selected physically-motivated criteria, such as a given parameter 
space of spectral-line width, depth and velocity. Further targeted 
study of H I absorption systems will allow us to improve these 
criteria. 
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5 SUMMARY 

We have conducted a search for 21 cm associated H I absorption 
towards 29 target sources from the AT20G catalogue (Murphy 
et al. 2010), within the redshift range 0.04 < z < 0.08 and 
at a velocity resolution of 7kms~ 1 . Since these sources are se- 
lected at high frequency, they are expected to be relatively young 
and represent the most recently-triggered radio AGN. We detect 
three associated absorption systems, of which two were previ- 
ously unknown, yielding a 10 per cent detection rate. This detec- 
tion rate is consistent with that seen for previous target sources 
at similar radio luminosities, and perhaps given the typical opti- 
cal depth sensitivity of our survey. However these interpretations 
should be treated with caution until we have obtained signifi- 
cantly larger sample sizes from future large-scale surveys (e.g. 
ASKAP-FLASH with a proposed ~ 150, 000 lines-of-sight). 

Of the three detected H I absorption lines, two exhibit sep- 
arate broad, shallow and narrow, deep components, while the 
third consists of a single broad, shallow component. These sepa- 
rate components are attributed to the relative orientation of both 
a galactic disc (low velocity gas) and a circumnuclear disc (high 
velocity gas) with respect to the line-of-sight (e.g. Morganti et al. 
2011). In future work we aim to follow up the 2 previously 
unknown detections with long baseline interferometric observa- 
tions, in order to spatially resolve the component features. 

Our search for associated H I absorption in this sample was 
aided by a method of detection using Bayesian inference, de- 
veloped by Allison et al. (2011) for the ASKAP FLASH sur- 
vey, to both detect absorption and obtain the best-fitting number 
of components of the spectral-line profile. While this method is 
a powerful tool for performing automated searches of spectral 
lines in low signal-to-noise data, we must be wary of the limita- 
tions produced by other signals that contaminate the continuum. 
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